Stochastic Path Integral Formulation of Full Counting Statistics 
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We derive a stochastic path integral representation of counting statistics in semi-classical systems. 
The formalism is introduced on the simple case of a single chaotic cavity with two quantum point 
contacts, and then further generalized to find the propagator for charge distributions with an ar- 
bitrary number of counting fields and generalized charges. The counting statistics is given by the 
saddle point approximation to the path integral, and fluctuations around the saddle point are sup- 
pressed in the semi-classical approximation. We use this approach to derive the current cumulants 
of a chaotic cavity in the hot-electron regime. 
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Noise properties of electrical conductors are interesting 
because they reveal additional information beyond linear 
response |l] . In the pioneering work of Levitov and Leso- 
vik ||] , the optics concept of full counting statistics (FCS) 
for photons was introduced for electrons in the context of 
mesoscopic physics. FCS gives the probability of count- 
ing a certain number of particles at a measurement ap- 
paratus in a certain amount of time and finds not only 
conductance and shot noise, but all higher current cumu- 
lants as well. Several methods have been used in finding 
this quantity. Originally, quantum mechanical methods 
based on scattering theory the Keldysh approach 

put forth by Nazarov [|| or sigma-models have been 
advanced and have been successfully applied to a number 
of problems among which we mention only multitermi- 
nal structures normal-superconducting samples 0, 
combined photon/electron statistics and conductors 
which are current (instead of voltage) biased |^ . 

A quantum mechanical treatment of transport shows 
that the leading contribution to current cumulants is of 
the order of the channel number N. For many conductors 
or circuits of interest, this leading order is a semi-classical 
quantity . Weak localization or universal conductance 
fluctuations provide only a small correction of order 1. 
Clearly, it is desirable to have a purely semi-classical the- 
ory to calculate semi-classical results. To provide such a 
derivation of FCS is the main purpose of this work. 

That a purely classical theory should be developed was 
realized by de Jong |ll] who put forth a discussion for 
problems which can be described with the help of mas- 
ter equations. A more general approach, leading to a 
set of rules for a cascade of higher order cumulants, was 
recently invented by Nagaev and applied to chaotic 
cavities p^ . The work presented here aims at provid- 
ing a foundation for the cascade approach by deriving a 
functional integral from which FCS, but also dynamical 
quantities such as correlation functions, can be obtained. 

The approach provided here applies to an arbitrary 
mesoscopic network. Its semi-classical nature does not 
allow the treatment of weak localization corrections nor 
is it applicable to macroscopic quantum effects like the 




FIG. f: Full counting statistics of a chaotic cavity at zero 
temperature: Comparison between hot (thick lines) and cold 
(thin lines) electron regime for both symmetric {{Q) = 40e) 
and asymmetric ((Q) = 50e) cases. Inset: The chaotic cavity 
with two point contacts. 

proximity effect. As in Boltzmann-Langevin theory [Q, 
our approach is based on a separation of time scales: fast 
microscopic fluctuations causes variations of conserved 
quantities (like the charge inside the conductor) on much 
longer time scales. Whenever such a separation of time 
scales is present in a stochastic problem, the formal- 
ism outlined here is useful. This includes for example 
normal-diffusive wires ]T2| , p^ , superconducting- normal 
structures outside the proximity regime 1 16 , as well as 
stochastic problems beyond mesoscopic physics . 

Path integral derivation. To introduce our path inte- 
gral formalism, we consider a simple example of electron 
transport through a chaotic cavity (see Fig. |l|) which is 
a large conductor connected to two metallic leads L and 
R through ideal point contacts. The leads are at dif- 
ferent chemical potentials fijj and fiL = fJ-R + ^l^, caus- 
ing a current / to flow through the cavity. The number 
of modes in each point contact iV^.i?, and the number 
of states eVup (np being a Fermi density of states in 
the cavity) are large, so that the transport is classical. 
We assume elastic transport (no energy relaxation), zero 
temperature, and chaotic electron dynamics in the cav- 
ity ||l8|. Then, the state of the cavity is described by 
only one variable fc which is independent of the elec- 
tron's coordinate, momentum and energy in the interval 
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fiR < e < @, 0- The average values of outgoing 
currents Il,r are given by (Ir) = {e^V/2'n:h)NR{fc) and 
{II) = — (e^F/27r?i)iVi(l — (/c)). Since currents are con- 
served {II) + (Ir) = 0, the average current through the 
cavity is (/) = {e''V/2TTn)[NLNR/{NL + Nr)], i.e. the 
resistances of the two point contacts add. 

The currents Il,r however fluctuate as a result of the 
discretness of the electron charge. These fast fluctua- 
tions (correlated on the time scale tq ^ h/eV) cause 
slow variations of the occupation fc on the scale of the 
relaxation time tc = 2TThnF/{NL + Nr) |2). The 
cause of this is charge conservation: II + Ir = —Qc, 
where Qc = e^npVfc is the charge accumulated in the 
cavity. Fluctuations of fc in turn affect the currents Il,r 
thereby generating correlations between them. Our goal 
is to integrate out currents Il,r taking into account cor- 
relations and to obtain the FCS of the transmitted charge 
Qit) - (1/2) J*dt'[lR{t') - hit')]. More precisely, we 
have to find the generator S{x, t) of the charge cumu- 
lants defined as 

P(Q,t) = {2Try' J dxeM-^xQ + Six.m (1) 
(Q'»(t)) = d"^S{x,t)/d{txr\x-o (2) 

where P{Q,t) is the distribution of transmitted charge, 
and {Q'^{t)) is its m-th cumulant. The average cur- 
rent and its zero- frequency noise power are given respec- 
tively by the first cumulant {Q{t))/t, and second cumu- 
lant {Q'^{t))/t, after taking the limit t — > oo. 

The separation of time scales tc/tq ^ 1 allows us 
to consider the intermediate time interval At, such that 
To <C At <C Tc, on which the occupation fc is approx- 
imately constant in time. Then the charges transmit- 
ted through the point contacts during this time interval, 
Ql,r = Jq dt'lL^R, are not correlated, i.e. the total dis- 
tribution is a product Pl{Ql)Pr{Qr)- On the other 
hand, since At is large compared to the correlation time 
To, their distributions Pl,r can be written as Fourier in- 
tegrals (|l|) with the cumulant generators proportional to 
time At |2^, namely Sl.r — Hh^R^t. They are well 
known |0] and for ideal point contacts given by 



Hl,r = {NL.,R/2TTh) j de{\n[l + fie){e''^-'^-l)] 

+ ln[l + /L.fl(e)(e-^^>^->«-l)]}. (3) 

For elastic transport, only the energy interval fiR < s < 
contributes with = 1, = 0, and f = fc- 
To proceed with our derivation we first discretize time, 
t„ = nAt, n = 1,...,A'^, and write the characteristic 
function of transmitted charge exp[S'(x,t)] for the un- 
conditional probability distribution 

exp[5(x,t)] = n n / dQlitn)Pimtn)] 
n l=L,R'' 

X exp(^^[Qfl(t„)-Qi(t„)]j. (4) 



Next, we impose a constraint 
fc{tn+i) = fcitn) - {e^nFV)-^[QL{tn) + Qi?(i„)], (5) 

which guarantees the conservation of charge on the cav- 
ity, and gives dynamics to the fluctuating charge Qc- 
Integrating out the charges QL,R{tn) in Eq. with the 
use of Eqs. (|l|) and (|^), and taking a continuum limit, we 
obtain one of our main results 



exp[^(x,t)] = (27r)-iy VQc^xc 



xexpjy^ dt'[ixcQc + H{x,xc,Qc)]j , (6) 

where H{x,Xc,Qc) is 

H = HlIxc ~ x/2, fc) + Hr{xc + x/2, fc) (7) 

with H^ R given by Eq. (|3|). Eq. (||) resembles the imag- 
inary time path integral for the evolution operator of a 
quantum particle with coordinate Qci momentum xCi 
and Hamiltonian —H with the difference that H given 
by (0) is not Hermitian. Nevertheless, we can apply the 
Hamiltonian formalism to our path integral Next, 
we derive the saddle-point solution and discuss its appli- 
cability. 

Saddle-point solution. The saddle point of the path 
integral (^) can be obtained by a variation of the action 
with respect to xc and Qc and gives "classical equations 
of motion" 



iQc 



= -dH/dxc, ixc = dH/dQc 



(8) 



These equations have to be solved with the conditions 
Qc = Qc{0) for t' = 0, and xc = for t' = t ||. The 
solution has to be substituted into Eq. (^ to obtain the 
saddle-point part 5*0 of the cumulant generator. 

Eqs. (||) describe the relaxation of the initial state 
Qc{0) to the stationary state {xc, Qc} given by the so- 
lution of dH/dxc = dH/dQc = 0. On time scales large 
compared to the relaxation time tc , the initial and flnal 
integration points contribute little (see discussion below) , 
and we can neglect the first term in the action (||). This 
gives us a large time asymptotics [p5|, 



So{x,t) ^tH{x,Xc,Qc), t:>Tc, 



(9) 



which will be analyzed below. 

To further simplify the analysis, we concentrate on the 
symmetric cavity with equal number of modes = 
Nr = Npc in the contacts. Then the stationary saddle- 
point solution can be found analytically giving xc — 0, 
fc = 1/2, and 

Soix, t) = {eVNp,t/7rh) ln[(l/2)(l + e'^^^^)] (10) 

which is the known result for a symmetric cavity found 
from quantum mechanical calculations p6|. 
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To investigate the validity of the saddle-point solutions 
and (H) we calculate the contribution of Gaussian fluc- 
tuations by expanding S in the vicinity of the stationary 
point to second order in xc ^^^d Qc ■ For the correction to 
So we obtain the action of a harmonic oscillator in imagi- 
nary time with inverse frequency w"^ ~ tc cos(ex/4) 
and mass M — (47rS/e^T^Afpc) cos^(ex/4). In the con- 
text of our problem such an action describes the linear 
dissipative dynamics Qc — ~^{Qc ~ Qc) + ^{t) with 
the Gaussian Langevin source {v{t)v{t')) — M~^6{t — t') 
p^t , which takes the initial state Qc(0) to the station- 
ary state Qc after time t ^ tc- Integrating out Qc we 
obtain for the total action S = Sq + SSq + Sfi, where 
SSq = — (l/2)Mti;(5p(0) tanh(wi) is the initial condition 
contribution and Sfi = — (1/2) ln[cosh(a't)] is the fluc- 
tuation contribution. Since in general Qc'(O) ^ e^npV, 
the part 5Sq is small compared to Sq ii t 3> tc- The 
part Sfi is small compared to So if enpV 3> 1, which 
has been assumed for our classical cavity All this 

justifies the saddle-point solution (||). Having established 
the framework of our formalism, we now proceed with its 
generalization and applications. 

Generalization. We consider the nonequilibrium dy- 
namics of an arbitrary classical system, which can be 
described by fast fluctuating currents Iaf3 and slow fluc- 
tuating charges Qa- Metallic reservoirs, such as leads 
and cavities, can be taken into account by replacing a 
(a, Z), where the index I — 1,...,L enumerates reser- 
voirs, and a enumerates generalized charges in the reser- 
voirs. These charges Qa are defined a,s Qa — ~ J2f3 ^ap- 
Assuming that the probability distributions of fluctuat- 
ing currents laji are known, we have to find an evo- 
lution of the distribution Y{Q,t) of the set of charges 
Q = {Qa} for a given initial condition T{Q^ 0). In other 
words, one has to find an evolution operator U{Q,Q',t) 
such that T{Q,t) ^ J dQ'U{Q,Q' ,t)r{Q' ,0). We use 
again the separation of time scales, follow the lines of 
the derivation of Eq. (^), and obtain the evolution oper- 
ator U — exp[So{Q,Q' ,t)], with Sq being a saddle-point 
solution of the action 



S^ f di'[zxQ+(l/2) ViJcMxc 

•^0 a/3 



(11) 



where x = {Xa} is the set of charge counting fields. The 
generators Hap of the cumulants of currents Iaf3 have the 
obvious symmetry HapiXa - Xfs) = H^aixfi - Xq)- 

Next we assume that the "Hamiltonian" H ~ 
{f/2)^^p Hap depends only on the subset of charges 
Q'^ = {Qa}, which we call conserved charges since they 
are related to conserved quantities, such as total energy 
and charge. We complete the set Q — {Q'^, Q°} by the 
subset of non-conserved or absorbed charges Q"" — {Qa}- 
Since H does not depend on Q", we find that the cor- 
responding "momentums" x° = {Xa} integrals of 
motion: x° = const. This allows us to integrate out the 



corresponding terms in ([ll|): i dt'x'^Q"' = ix°[Q°(0 — 
(5°(0)]. Thus the rest of the action (|l|) becomes a gener- 
ator of cumulants of the absorbed charge Q°(t) — Q"(0). 
In particular, the stationary solution for the cumulant 
generator is 5*0 = tH {x"" , x'^ , Q'^) , where x'^ and Q'^ are 
solutions of the equations dH/dx'^ — dH/dQ'^ — 0. 

In the example of the elastic transport considered 
above there are 3 reservoirs: Two leads {a — 1,3) and 
one cavity (a = 2). The state of the cavity is described 
by only one charge Qc — Q2 , the total charge on the cav- 
ity, while Q = (Qa — Qi)/2 is the absorbed charge which 
is being counted. The counting fields are xc = X2, and 
x/2 = Xi = ~X3- Next we consider the case of inelas- 
tic transport, which requires the introduction of a new 
generalized conserved charge Ec , the total energy of the 
cavity. 

Hot- electron regime. We consider the electron trans- 
port through a cavity in the so-called hot-electron regime, 
i.e. when the electron-electron scattering time is much 
shorter than the relaxation time Tee <SC tc, electrons in 
the cavity then relax to local thermal equilibrium before 
escaping from the cavity. Their distribution is given by a 
Fermi function /(e) = {l+exp[{e—iJc)/Tc]}~^ where /ic 
and Tc denote respectively the electro-chemical potential 
and electron temperature in the cavity. The noise in the 
hot-electron regime was first considered theoretically for 
diffusive wir es p9| . This calculation and its experimental 
verification ||3C| found that electron-electron interaction 
increases shot noise. Similar results for the chaotic cav- 
ity were presented in Ref . ||T^ and measured in Ref . . 
We now derive the FCS of transmitted charge Q. 

Our calculation starts from the observation that the 
nonequilibrium state of the cavity is fully described by 
only two parameters, nc and Tc. They can be ex- 
pressed in terms of conserved values: The total charge 
Qc — e^^Cf^fic and the total energy Ec = {l/2C'^)Q'^ + 
QcVg + {t^"^ /^)nFTc of the cavity, where C^j, is the elec- 
trochemical capacitance of the cavity, 1/C^ = 1/C -\- 
l/{e^np), and Vc is the gate voltage. The charge and 
energy conservation can be written as 



de[lLis) + lRie)]^~Qc, (12) 
dee[lL{e) + lR{e)]^-eEc, (13) 



where Il.r.{£) are the outgoing currents through the 
point contacts per energy interval de. We now replace 
the constraint (^ with these two equations, and apply 
it to the energy resolved analogue of path integral (|). 
By doing so, we introduce the counting field xe for the 
total energy Ec. Instead of (0), the new "Hamiltonian" 
H{x,Xc,XE;fJ.c,Tc) takes the form 

H = J de[HL{xc + £XE-x/2) + HR{xc + exE + x/'2)], 

(14) 
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with Hl,r given by Eq. (||). 

To obtain the FCS in the long time Hmit t ^ tc, where 
Tc = 2TinC^/[e^{NL + Nr)] is the RC-time (the time of 
charge screening in the cavity |22j), we need to evalu- 
ate H as Si function of x at the saddle point. We find 
the saddle-point solution S^){x^t) numerically and calcu- 
late P{Q,t) by Fourier transformation (|l|). The result is 
shown in Fig. |^. It is clearly visible that noise is enhanced 
in the hot-electron regime. 

In order to make further analytical progress we note 
that the saddle point solution for x = is given by the 
chemical potential ^(7 = (^LAiL + ^7?,Mfl)/(-^L+^_R) and 
the effective temperature Tc = \/T'^ + Ty , where T is 
the temperature of the reservoirs, and Ty — {^/i/'K)rieV 
is proportional to bias with -q^ = NlNr/ {Nl + NrY. 
The internal fields XCtXe are zero for x = 0. It is now 
straightforward to expand Eq. (|lj) and the correspond- 
ing saddle point equations around x = and to calcu- 
late higher cumulants order by order. This procedure is 
very similar to the cascade approach considered in Refs. 
WA &3|. For the first few cumulants we obtain 



2TTh 



VNlNr (T + Tc) ., 



ttTc 

T-Tc 



2T^ 



(15a) 
(15b) 
, (15c) 



where the second cumulant coincides with the one in the 
Ref. The result for the third cumulant shows that 

odd cumulants are strongly suppressed. We note that 
screening does not affect the FCS in the long time limit 
t > rc. 

Conclusions. We have constructed a stochastic path 
integral formulation of full counting statistics. Our ap- 
proach is based on the fact that fast microscopic quan- 
tum fluctuations give rise to slow variations of conserved 
quantities. The method has been illustrated with chaotic 
cavities and novel results have been presented for the 
hot-electron regime. We emphasize the general nature of 
this approach and its applicability to stochastic problems 
even outside mesoscopic physics. 
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National Science Foundation. 
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